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Abstract 



Using a quadratic saddle-point approximation, we show how information 
about a particle-emitting source can be extracted from gaussian fits to two- 



particle correlation data. Although the formalism is completely general, ex- 



^ ■ traction of the relevant parameters is much simpler for sources within an 



interesting class of azimuthally symmetric models. After discussing the stan- 
dard fitting procedure, we introduce a new gaussian fitting procedure which 
is an azimuthally symmetric generalization of the Yano-Koonin formalism for 
spherically symmetric sources. This new fitting procedure has the advantage 
that in addition to being able to measure source parameters in a fixed frame 
or the longitudinally co-moving system, it can also measure these parameters 

in the local rest frame of the source. 
PACS: 25.70.Pq 
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I. INTRODUCTION 



Recently quite a bit of work has been done in trying to determine which attributes of 
the hadronic source formed in high-energy particle or heavy-ion collisions can be determined 
by measuring the Hanbury-Brown-Twiss (HBT) correlations of identical emitted particles. 
Usually the experimental correlation function is fit with a gaussian in some components 
of the four-momentum difference q = pi — P2 (^Hil- The parameters of such a fit (called 
correlation radii) are then often compared to some simple analytic model in order to get an 
idea of what is being measured. For example, for a static gaussian source, the duration of 
emission time is directly proportional to the difference of the squares of the correlation radii 
which are parallel to ("out") and perpendicular to ("side") the transverse component of 
the total pair momentum. However, this result is not true if for example the actual source 
contains any z-t, z-x, or x-t correlations, such as would be caused for quickly expanding 
sources. 

The purpose of this paper is to determine exactly what features of the source are actually 
being measured by experimental correlation data. To do this, we use a quadratic saddle- 
point approximation |§-|8| to a general source function in order to derive a completely Lorentz 
covariant expression for the two-particle correlation function which can be applied to a wide 
range of analytic models. In this approximation, 10 K-dependent parameters are needed 
to describe a general source, where K = |(pi + P2) is the average momentum of the two 
particles. The origin of these 10 parameters can be understood by noting that for each value 
of K, our approximation is mathematically equivalent to a gaussian ellipsoid described by 
three Euler angles of orientation, three components of the velocity of the local rest frame, 
three spatial extensions, and one temporal extension. As we will show, however, only six 
K-dependent parameters (not including the chaoticity parameter A) can be measured by 
making a gaussian fit in q to the correlation function. 

Furthermore, a source function which is azimuthally symmetric in coordinate space will 
not in general be azimuthally symmetric in momentum space, since the direction defined by 



2 



K± breaks this symmetry. As a result, the K-dependent ellipsoid equivalent to such a source 
still requires one Euler angle, two velocity components, three spatial extensions, and one 
temporal extension in order to describe it. In this case, only four K-dependent parameters 
can be determined by making a gaussian fit in q. Since for arbitrary K± the number of 
source parameters exceeds the number of gaussian-fit parameters, some definite model must 
be used in order to interpret how the latter depend on the former. 

For pairs with K± = 0, however, only four K-dependent parameters are needed to 
describe the source (one velocity component, two spatial extensions and one temporal ex- 
tension), so the four parameters measured in a gaussian fit provide enough information to 
unambiguously determine all of these source parameters. Although in practice it is very 
difficult to measure a pair whose total transverse momentum vanishes, there is an interest- 
ing class of models in which the K± = simplifications persist for nonzero values of K±. 
In Sec. 5 we study a hydrodynamical model similar to ones which have been used to fit 
one-particle distributions from heavy-ion collisions at the Brookhaven Alternating Gradient 
Synchrotron (AGS) and the CERN Super Proton Synchrotron (SPS). Within the context 
of this model, one would expect the simplifications in the extraction of source parameters 
to occur for pairs whose average energy in the measurement frame is less than the freezeout 
temperature divided by the square of the transverse flow velocity (about 560 MeV for Si+Au 
at the AGS or S+Pb at the SPS). 

If the correlation function for a model in the afore-mentioned class is fit to a gaussian 
in the spatial components of q, then in order to extract the relevant source parameters it is 
crucial to pick beforehand the correct longitudinal reference frame for the measurement. If, 
on the other hand, a fit is made to a generalization of the Yano-Koonin formalism p|,|T0" 



the correlation radii will automatically measure the relevant source parameters even if the 
wrong frame is chosen for measurement. For example, for a finite longitudinally expanding 
source, the new formalism allows measurement of the source parameters in the local rest 
frame of the fluid, whereas fixed and longitudinally co-moving system (LCMS) radii measure 
the source in different frames. 



II. LORENTZ CO VARIANT CORRELATION FUNCTION 



For the two-particle correlation function, we use the well-established theoretical approx- 
imation 

C(a K)~l± ljggjfeg^ (21) 

where q = Pi - P2, % = E x - E 2 , K = |(pi + p 2 ), and K = E K = ym 2 + |K| 2 . The 
plus sign is to be used for boson pairs and the minus sign for fermion pairs. The labeling 
of particles 1 and 2 is defined such that q\ = q x is always positive. In this way, pairs with 
positive q2 = q y and/or q 3 = q z are physically distinct from those with negative q 2 and/or 
q 3 . The S(x, K) in Eq. ( |2.1|) is a function which describes the phase-space density of the 
emitting source. 

The spacetime saddle point x(K) of the emission function S(x, K) is defined via the four 
equations HH 





\nS(x, K) _ = , (2.2) 



Ox,, 

where /x = {0, 1, 2, 3}. The saddle point is that point in spacetime which has the maximum 
probability of emitting a particle with momentum K. A quadratic saddle-point approxima- 
tion for S(x, K) then yields 

S(x,K) ~ S(x,K)exp [-\{x - xY{x - x) u B^{K)\ , (2.3) 

where the symmetric curvature tensor B^ v is given by 

B^iK) = -d^d u \nS(x, K) \_ . (2.4) 

We define the curvature radius in the /ith direction by ||[7j 

A„(K) = [B„(K)}- 1/2 . (2.5) 

Note that since B^ V {K) is symmetric, it will in general have 10 independent components. 
From the form of Eq. (|2.3|) it is seen that the saddle-point approximation is mathematically 



equivalent to an ellipsoid described by the 10 K-dependent parameters mentioned in the 
Introduction. As long as the saddle point x(K) is unique, knowledge of the 10 functions 
B^ V {K) is in most practical situations sufficient for a complete characterization of the source. 
It is convenient to define p,[L3],[14| the following K-dependent average of an arbitrary 



spacetime function with the source density S(x, K): 



Using this notation, the correlation function (|2.1|) can be compactly written as 



C(q,K) ~ 1 ± (e iq - x )(K) . (2.7) 
Furthermore, within the saddle-point approximation ( |2.3| ), the following relations hold: 

W = x,(K) , (2.8) 

{x„x v ) - (x M > (x v ) = (£~V( K ) . (2.9) 

The saddle point is thus the average spacetime point from which particle pairs with momen- 
tum K are emitted, and the components of the inverse of the curvature tensor (S _1 ) M!/ (K) 
give the spacetime correlations of the source. The four diagonal elements 

(S-^WCK) = (xl) - (x,) 2 (2.10) 

can be understood as the squares of the lengths of homogeneity of the source as seen by 
pairs with momentum K. It should be noted that the homogeneity lengths agree with 
the curvature radii (|2.5| ) only if the curvature tensor B^ v is diagonal. This was implicitly 
assumed by the authors of ||[7|], who first introduced the name "homogeneity length" but 
used it for the curvature radii ( |2.5|) of the source near the saddle point. 

Within the approximations of Eqs. (|2.1| ) and ( |2.3|) , calculation of a general correlation 
function is straightforward, yielding 

C(q,K) = 1 ±exp l-q^iB- 1 ) J . (2.11) 



From Eq. (|2.9|) we can see that the correlation function directly measures the spacetime 
correlations within the source. 

It may at first seem that all of the components of the correlation tensor B~ x can be 
found simply by comparing the results of a four-dimensional fit to the correlation function 
with Eq. fl2.11| ). Such a fit is not possible, however, since q is highly correlated with the 
other components of q through the equation 

2K, 



ft) = X) Aft j where # 



E l + E 2 



[2.12) 



By making the approximation 



A ^ Ki/E K 



(2.13) 



which is valid for pairs with |q| <C Ek, one can use ( |2.12 ) to fit the correlation function 
( [TTI| ) to the form 



C(q,K) = 1 ± Aexp [-qlR\{K) - q 2 2 R 2 2 (K) - q 2 3 R 2 3 (K) 

-2q iq2 R 2 12 (K) - 2 qi q 3 R 2 13 (K) - 2q 2 q 3 R 2 23 (K) 



(2.14) 



where the R 2 j cross terms can be either positive or negative and A is a parameter introduced 



to allow for coherence effects []15| , [l6|| and/or particles from the decay of long-lived resonances 
17|-p0|. From Q2.9|) and ( [2.1 2|) the six functions i??(K) (Ri = Ru) can be expressed as the 



correlations |qjl4| 



R 2 AK) = ((xi- Pit)(xj - pjt)) - ( Xi - ^t)(xj - /3jt) . 



(2.15) 



In general, the six R 2 parameters found by fitting correlations to Eq. (|2.14j) do not provide 



enough information to determine the 10 independent components of (B 



Due to Eqs. ( |2.12 ) and (|2.13 ), any Lorentz transformation and/or spatial rotation to a 
new coordinate system can be written as a purely spatial linear transformation, 



Qi = aijqj 



(2.16) 
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Thus the new R' 2 parameters found in the primed frame will simply be linear combinations 
of the R 2 parameters found in the original frame. For example, the longitudinally co-moving 
system (LCMS) is defined as the longitudinally boosted frame in which j3' 3 = P,|5|,pT|| . The 
R' 2 parameters in this (primed) frame are related to those in some fixed (unprimed) frame 



via 



R' 2 = R\ + ^PlPlRi + 2 7 2 /3 1 /? 3 i??3 
R' 2 2 = R\ + i A 0\0lR\ + 2 7 2 /? 2 /? 3J R 2 2 3 
R'i = 7 2 i? 3 2 

R'12 = RI2 + i 2 mR 2 13 + rVW^s + I'PiWlRl 

R' 2 \ = lRlz + l*(32foRl, (2.17) 



where the /3j are measured in the fixed frame and 7 = 1/y 1 — [3 2 . The above equalities can be 
used as an experimental test of the validity of the saddle point approximation in the following 
way: If measured LCMS radii are not equal to the above corresponding combinations of 
fixed-frame radii, then the saddle point formalism is not a good approximation to the actual 
correlation function. 

For sources which have highly non-gaussian spacetime dependencies, it is better to define 



[22 1 the inverse of the curvature tensor directly through Eqns. ( |2.9| ) and ( |2.6| ) (using the full 
source S(x,K)) rather than to use second derivatives to define the curvature tensor at the 
saddle point via fl2.4p . A simple example of when this is necessary is provided by the emission 



function from a uniform sphere: 

S(x, K) = f(K)5(t - t )6(R - r) , (2.18) 

where r = y/x 2 + y 2 + z 2 and R is the radius of the sphere. For this source function, the 
saddle point is not unique and second derivatives are ill-defined. Nevertheless, (B^ 1 )^ in 
( O ) is perfectly well defined, producing through ( |2.11| ) a gaussian correlation function of 



the form 



C(q, K) = 1 ± exp (-|q| 2 i? 2 /5) • 



(2.19) 



The reader can verify that the above expression is a good approximation to the correlation 
function 
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C(q,K) = 1 ±9 



cos(|q|i?) sin(|q|i2) 



(2.20) 



(HR) 2 (HR) 3 

which is found by plugging ( 2.18]) directly into (|2.1|) . This example clearly shows that in 
terms of the correlation function, even a highly non-gaussian source can be well approximated 
by a gaussian with the same rms width. 



III. AZIMUTHALLY SYMMETRIC SOURCES 

For an azimuthally symmetric source, it is convenient to choose z to point along the 
beam ( "longitudinal" ) axis and to choose x to point in the same direction as the component 
of K which is perpendicular to the beam ("out"). The remaining ("side") direction is then 
defined by y = z x x ||23|| . Note that by definition K\ and (3\ are always positive and 



K 2 — fa — 0. Since the latter is true, azimuthally symmetric sources must satisfy 

S(t, x, y, z, K) = S(t, x, -y, z, K) . (3.1) 

From Eq. ( |2.3|) this implies that y = B^ 2 — B 2f _ l = for /i 7^ 2, so for azimuthally symmetric 
sources B^ u has only seven independent components. 

By inserting Eq. ( |3.1| ) into Eq. ( |2.1| ), one can see that the correlation function from an 
azimuthally symmetric source must be unchanged under the substitution q 2 — > —q 2 - This 
implies that R\ 2 = R 23 = 0, so that only four R 2 parameters can be found by making a 
gaussian fit to the correlation function [HIT 



C(q,K) = 1 ± Xexp[-qiRi(K) - q 2 2 R 2 2 (K) - q 2 3 Ri(K) - 2 qi q 3 Rf 3 (K)\ . (3.2) 

Expressing these four correlation radii in terms of the seven independent elements of the 
curvature tensor we have 
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Rl = r A 2 (l - A^i^) + (3f\l(l - XlXlBl) + 2ftA?Ajj(B 10 - A^ 3lJ B 30 ) 



-^3 = r Ag(l - AqA^Si ) + Pi Aq(1 - A^AgS^) + 2/? 3 A 3 Ao(I? 3 o - A^S 31 S 10 ) 



-^13 = r /3i/?3Ao(l - A^AgS^) + /^AgA^-E^o - AiS 31 S 10 ) 



+/?3AiAo(Sio - A 3 £? 3 i-B 30 ) - \\\\(B2,i - XlB 10 B 30 ) 



(3-3) 



where 



F— 1 \ 2 \ 2 R 2 \2\2n2 \2\2 R 2 , o\2\2\2d R R 

1 - i - A A 1 IS 10 - A X A 3 £S 31 - A A 3 # 30 + ZA A 1 A 3 &i tf3otf: 



(3.4) 



It is now apparent that the statement that R\ — R\ is proportional to the square of the 
emission time is highly model dependent. Implicit in this statement are the assumptions 
that Ai = A2 and that all t, and z-t correlations are negligible. 

Until now we have been discussing fits to the correlation function which are three- 
dimensional gaussians in the spatial components of q. It is also possible to fit correlation 



functions in a different three-dimensional space defined by q± = Jq\ + g§, ^3 and go- As 
we pointed out in the last section, go is highly correlated with the spatial components of q 
through Eq. (|2.12|) which we rewrite here as 



where is the angle between x and q^. Due to this correlation, g only varies over a finite 
range which is nonetheless greater than zero for any nonzero values of Pi and q±_: 



Despite this phase-space limitation, it still possible to fit correlation functions with gaussians 
in a (g±,g 3 ,go) space JZ5|]. 

For example, one could use an azimuthally symmetric generalization of the Yano-Koonin 
formalism [3, |10H 




Qo = (3\q± cos + (3 3 q 3 , 



(3.5) 



- Piq± + p3Q3 < go < Piq± + PsQ3 ■ 



(3.6) 



C(q, K) = 1 ± A exp -q\R\ + (g 2 - g 3 2 )i? 2 - (g-f/) 2 (i? 2 + R\) 



(3.7) 



where 

I7 = 7(l, 0, 0, v) , 7= 1/Vl-v 2 (3.8) 

(in units with c = 1), and i? t (K), i? (K), i? 4 (K), and i>(K) are the four fit parameters. When 
using this fitting procedure, it is convenient to define the particle labeling such that go is 
always positive. In this way, pairs with positive and negative ^3 are physically distinct and 
can be separately binned. For the remainder of this paper, we will refer to Eqs. ( |3.2|) and ( |3.7|) 
as the "standard" and "GYK" (generalized Yano-Koonin) fitting procedures, respectively. 

The boost-invariant form of the "GYK" fit means that for any given value of K there 
is a longitudinally boosted reference frame in which u(K) = 0. The fit parameters Rt(K), 
Rq(K), and i? 4 (K) measure the source in this frame, regardless of which longitudinally 
boosted frame is chosen for the evaluation of the g M . Although the general interpretation of 
the v(K) = frame is difficult, it can be shown that in this frame 

((Piz - MWit - x)) - (fcz - fax)(frt -x)- (3 3 (y 2 ) = . (3.9) 

Furthermore, in this frame the remaining fit parameters take the form (for f3\ 7^ 0) 



Rl = (y 2 ) = A 



Pa_\ \ I Pa_\ &/ 2 



Rt = ( Z-^x - (z-^x) -%(y 



-'\ / , \ 2 ^2 

,1 a. 1 V\ A 1 \ 2 I/..2 



R °=\{ t -r)/-\ t -r) (3 - 10) 

Just as for the correlation radii of ( |2.15|) and ( p.3| ), these fit parameters contain mixtures 



of various spatial and temporal lengths of homogeneity, making the extraction of source 
parameters in general highly model dependent. However, a "GYK" fit has the advantage 
that the time structure of the source enters in only one of the radius parameters, namely 
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IV. THE MEASUREMENT FRAME 



In order to illustrate the utility of the saddle-point formalism and the importance of 
picking the right longitudinal reference frame for making measurements, we will now study 
three example source models. 

A. Static Gaussian Source 

Consider a source function which is defined by spacetime gaussians in its center-of-mass 
(CM) frame: 

x 2 + y 2 z 2 it -t a ) 2 ' 



S(x J K) = f(K)exp 



(4.1) 



2R 2 2L 2 2(At) 2 

Using Eq. ( |2.2| ), we see that the saddle point for this function is independent of K and given 
bjx = y = z = and i = to. In the CM frame, all of the off-diagonal components of 
vanish, while the curvature radii are given by 

\ 1 = \ 2 = R , A 3 = L , and X = At . (4.2) 

That the curvature radii are independent of K is true for any source in which the spacetime 
and momentum dependences factorize. For the simple source of Eq. fl4.1|) they simply 
measure the relevant geometrical "radii" of the system. 

If one makes a "standard" fit to the correlation function using CM momentum differences, 
the extracted radii will have the interpretation 

R 2 1 = R 2 + (3 2 (At) 2 , R\ = R 2 

R\ = L 2 + Pl(At) 2 , R 2 n = (3MAt) 2 (4.3) 

The duration of emission (At) can thus be extracted either from R 2 3 or from the difference 
R\ — R\. If some frame other than the CM is used for making a "standard" fit, then the 
correlation radii will have much more complicated dependencies on the source parameters. 
In particular, R\ — R\ will not be proportional to (At) 2 . 
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If a "GYK" fit is made, on the other hand, then in the CM frame or any frame longitu- 
dinally boosted from it, one will find the same expressions for the correlation radii: 

R t = R, R4 — L , R = At . (4.4) 

Furthermore, the parameter v will measure the velocity difference between the CM and 
measurement frames. In other words, even if measurements are made in the "wrong" frame, 
a "GYK" fit will still produce correlation radii which measure the source in its rest frame. 



B. Boost-Invariant Source 



Next we consider a longitudinally expanding, boost-invariant source: 



S{x,K) 



(At) 



cxp 



K-u x 2 + y 2 (t - r ) 



T 



2R 2 



2(Ar) 2 



(4.5) 



Here T is a constant freeze-out temperature, r = \/t 2 — z 2 is the longitudinal proper time, 
and m t = \jm 2 + Kj_. The longitudinal expansion of the source is described by the flow 
four-velocity relative to some fixed frame 



■u M = ( ch?7, 0, 0, shi] ) , 



(4.6) 



where t\ — \ ln[(t + z)/(t — z)\ is spacetime rapidity. 

Calculation of the saddle point in the fixed frame is straightforward, yielding x = y = 0, 
t = t and fj = Y, where Y is the rapidity of a particle with momentum K. Unlike in the 
previous model, both the saddle point and the curvature tensor depend on K. Explicitly, 
the curvature radii and off-diagonal elements of B^ u are given by 



A 
Ai = A2 

A 3 



sh 2 Y /m t \ ch 2 Y 



R 



~ch 2 Y 




sh 2 F " 


T 2 


V T j 


+ (Ar) 2 _ 



-1/2 



-1/2 



5, 



03 



B 30 = -shY chY 
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(Ary 



(4.7) 



Here we see that the z-t correlations of the longitudinally expanding source give rise to a 
nonvanishing B 30 for pairs with 

The correlation radii from a "standard" fit can be found by plugging these source pa- 
rameters into Eq. (|3.3|). One finds 



p2 p2 

1X2 — -ft 



R\ = R 2 + (3f 



(Ar) 2 ch 2 F + (J^) r 2 sh 2 y 



m t J ch 2 F 

R 2 13 = -(3i(— )r 2 thF. (4.8) 

Using a saddle point approximation to a hydrodynamic model similar to the one presented 
here, Makhlin and Sinyukov first derived the above expression for R\ [ 2E ] . In a recent paper 



by the NA35 collaboration at CERN, the mt and Y dependence of R3 measured in relativistic 
heavy-ion collisions was compared to that expression in order to estimate a freezeout proper 
time t [[|. This treatment was not consistent, however, since the parameter i? 2 3 was omitted 
from the fits. 

Notice that R\ — R\ in Eqs. ( |4.8|) is in general proportional neither to (At) 2 nor to A 2 ,. 
This is simply a result of choosing the wrong measurement frame. The extraction of At 
from correlation radii is greatly simplified if instead of measuring the correlation in a fixed 
frame, one measures it in the local rest frame of the source near its saddle point x(K). At 
that point, the flow velocity relative to the fixed frame is given by 

u = y chfj, 0, 0, shfj ^ (4.9) 

Since fj — Y, transforming to the local rest frame (primed) of the saddle point can be done 
for each pair by making the following substitutions 

rj' — rj — fj — rj — Y, Y' — Y — fj — 



q' 3 = chY(q 3 - thr g ) , % = chY(q - thY q 3 ) . (4. 10) 
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In other words, the local rest frame of a boost-invariant source is just the LCMS (Y 7 = 
frame) . 

It is easy to verify that in this frame B' 30 vanishes and 

\' = At, \' 1 = \' 2 = R, \' 3 = T ^T/m t . (4.11) 

Due to the diagonal nature of the curvature tensor, these source parameters can be easily 
extracted from the LCMS correlation radii 

Rf = R? + ^(Ar) 2 , R' 2 2 = R 2 , R' 3 2 = (T/m t )r 2 , Rf 3 = . (4.12) 

Note that in this frame R\ — R\ is proportional to (At) 2 . Using the "standard" fitting 
procedure, it obviously makes the most sense to measure correlations for boost-invariant 
sources in the LCMS rather than a fixed frame. 

Using the "GYK" fitting procedure, on the other hand, one finds that regardless of the 
longitudinal frame chosen, 

R t = R, R A = Toy/T/mt , R = At , (4.13) 

while v measures the difference between the measurement frame and the local rest frame 
of the source at the saddle point. For example, a measurement made in the LCMS would 
yield v = 0, while one made in a fixed frame would yield v = thY for pairs with an average 
rapidity Y relative to that frame. 

C. Finite Expanding Source 

The main problem with a boost-invariant source is that it gives rise to a dN/dy which 
is completely flat, whereas the dN/dy for produced particles which are actually observed in 
relativistic heavy-ion collisions are much better described by gaussians in rapidity. Further- 
more, in a boost-invariant source, Rf 3 will vanish in the LCMS, but nonzero values for Rf 3 
have been measured by NA35 |27| . These inconsistencies with boost-invariant sources lead 
us to consider an expanding model with a non-boost-invariant cutoff in spacetime rapidity 
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f8 yi^j28|j29|1 . Using the expansion four- velocity of Eq. ( (4.6|) , we define a source in its CM 
frame by 



(27r)y27r(Ar) 2 



K-u x 2 + y 2 (r — r f 



T 



2R 2 



(4.14) 



2(Ar) 2 2(Ar]) 2 _ 

In || it was shown that the At] cutoff term leads to more realistic gaussian-like rapidity 
distributions. The prefactor was introduced so that in the limit as At — > 0, S(x, K) becomes 
the Boltzmann approximation to a hydrodynamical source which freezes out at a constant 
temperature T and proper time r |30] . 

The saddle point for this model still has the coordinates x = y = and f = r , but the 
spacetime-rapidity coordinate for a given value of K (relative to the CM frame) is now given 
by the solution of 



. 



(4.15) 



It is apparent that the saddle point is only located at fj — Y for an infinite longitudinal tube 
At] —>■ oo or for pairs with Y = 0. 

As in the previous model, transforming to the local rest frame (primed) at the saddle 
point can be achieved by 



Tj —7] — Tj 



Y' = Y — fj 



<?3 = chr?(g 3 - thrjqo) 



q' = ch.7](q - thr]q 3 ) . 



(4.16) 



Notice that since Y' ^ 0, the local rest frame does not coincide with the LCMS. It is possible 
to show that in the local rest frame the curvature tensor becomes diagonal (B' 3Q = 0) and 



A' 3 = t 



\[ = \' 2 = R , X' Q = At 



T Ch ^- Y) - c h 2 (^Y) + lA^) 2 \ 



-1/2 



(4.17) 



In the CM and LCMS frames, however, the curvature tensor is in general not diagonal and 
each component becomes much more complicated. These complications carry over into the 
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correlation radii. For example, R\% does not vanish in these frames, and R\ — R\ is not 
proportional to either (At) 2 or Aq PJI4"||. 

By looking at one-particle slopes and rapidity distributions, it may be possible to estimate 
T and Arj. Using these values, one can numerically solve Eq. ( |4.15| ) for fj(K), use Eq. ( [4.1 6| ) 



to transform to the local rest frame, and then extract the source parameters by making a 
"standard" fit to the correlation function. Alternatively, one could simply make a "GYK" 
fit in the CM frame (or any other longitudinally boosted frame), and the result would be 

Rt = A' : = A' 2 , i?4 = A3 , Ro = Aq . (4-18) 

Again the parameter v (K) would measure the velocity difference between the measurement 
frame and the local rest frame. 

Each of the three models discussed above exhibits a different local rest frame. If one 
makes a "standard" fit to the correlation function, it is important to guess the correct refer- 
ence frame before performing the fit in order to extract useful information about the source. 
This implies that one must have some a priori knowledge about the source before making the 
fit. If one makes a "GYK" fit, however, there is no need to pick a frame beforehand. In all 
three of the above cases, the fitting procedure itself automatically chooses the correct frame 
in which to measure the source. This feature of the "GYK" fitting procedure is actually 
common to a whole class of models we will discuss below. 

D. A Class of Models 

All of the models we have discussed so far belong to a class of source functions which 
satisfy the following conditions: 

Ai = A 2 and B w = B 31 = . (4.19) 

For models in this class, the curvature tensor (|2.4|) takes the simple block-diagonal form 
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B 



[Ms 



Xn 2 



A 2 2 



V 



-B30 A 2 



(4.20) 



A 3 - z B 30 


Moreover its inverse, the correlation matrix B~ x of Eqs. ( [2.1 1| ) and ( |2.9| ), has the same 
block-diagonal structure, i.e., the x-t and x-z correlations (S _1 )io and (_B _1 )i 3 vanish, and 
(B- 1 ^ = (B~% 2 . 

Since models in this class have only four non-vanishing components of the curvature 
tensor, these four can be unambiguously determined by measuring the four parameters 
coming from either a "standard" or "GYK" fit. Explicitly, 
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^(R 2 + v 2 Rl) 

Z?2 p2 

rt /t 4 

^(Rj + v 2 R 2 ) 
V1 2 {R 2 + R 2 ) 

f?2 p2 



(4.21) 



where 7 is defined in (|3.8|). If the q used for fitting the correlation function are evaluated 
in a fixed frame, then the above parameters describe the source as seen in that fixed frame. 
If, on the other hand, LCMS values for are used, then Eqs. (|4.21|) (with f3 3 = 0) will 
determine LCMS source parameters. 

An interesting experimental test is provided by comparing the left and right sides of 
Eqs. ( |4.21| ). Namely, if after making both "standard" and "GYK" fits to the correlation 
function it is found that the left and right sides of Eqs. (|4 . 2 1|) are not equal, then the source 



in question cannot belong to the class ( |4.19| ). 

As we mentioned previously, the form of the "GYK" fit means that the extracted pa- 
rameters -Rj(K), i?4(K), and Rq(K) naturally measure the source in the v (K) = frame. 
Setting v = in Eqs. ( ^.21| ), we see that the right halves of the first three equations reduce 
to Eqs. ( |4.18| ). In other words, the R parameters directly measure the curvature radii of 
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the source in the v(K) = frame. Furthermore, the last equation in ( |4.21| ) shows us that 
the v (K) = frame corresponds to the frame in which B 30 (K) vanishes (z-t correlations 
vanish) and the curvature tensor diagonalizes. For many interesting models including the 
three discussed above, this frame also corresponds to the local rest frame of the source at 
the saddle point x(K). 

As for the "standard" fitting procedure, the left halves of Eqs. (|4.21| ) tell us that only 
if we are clever enough to pick the -B30 = frame beforehand will Aq reduce to the more 
familiar form 

\1 = ^ 2 {R\-Rl). (4.22) 

Pi 

Using RQMD events in the LCMS, this equation was recently shown to be a good approxi- 
mation only for pairs with very small Y and K± in the center-of-mass frame of a symmetric 



projectile-target collision |fL9|| . The restriction of small Y found by these authors can be 
explained in the following way: Only at Y = do the LCMS and the local rest frame defi- 
nitely coincide for all sources resulting from a symmetric projectile-target collision. As \Y\ 
increases, these two frames may begin to diverge just as they did for the model discussed 
is subsection |LV CJ above. For large Y, B 30 becomes non-negligible in the LCMS and ( f4.21| ) 



rather than (f4.22|) must be used to extract Aq. As we will show in the next section, the 



additional restriction of small K± can be explained if a source is undergoing transverse as 
well as longitudinal expansion. 

The reader should note that there is a difficulty with the "standard" fitting procedure 
which arises due to the fact that experimental correlation functions are always generated for 
pairs which lie in certain bins in K± and K 3 =Kl rather than for exact values of these average 
momenta. In extracting Ao and B 3Q of Eqs. ( |4.21| ) from the left-hand expressions, there will 
be a certain ambiguity as to which values of p\ and j3 3 should be used. Fortunately, this 
problem can be circumvented by employing the new fitting procedures that we introduce in 
Appendix A. Since no /3j are present on the right side of Eqs. ( (4.21[ ), the difficulty does not 
arise at all when one uses the "GYK" fitting procedure. 



Using the generalized Yano-Koonin formalism has the further advantage that by making 
a single fit in a fixed frame, one can determine the source parameters both in that frame via 
( [4.21| ) and simultaneously in the v (K) = frame via ( |4.18| ). Similarly, by making a single 
fit in the LCMS frame, one can simultaneously determine both LCMS and v (K) = source 
parameters. Using the "standard" fitting procedure, on the other hand, requires at least 
two fits (see Appendix A) just to determine the source parameters in a single frame. 

Finally, we would like to discuss the reason why we chose to study the class of models 
defined by Eqs. ( j4.19|) . First of all, as we have seen explicitly, there are a number of 
interesting models which naturally fall into this class. More importantly, however, the 
source function for pairs with K± = from any azimuthally symmetric model will always 
fall into this class. The reason for this is simply because when K± = 0, there is no way 
to distinguish between the "side" and "out" directions. Consequently, Ai = A 2 and since 
B211 = 0, it must also be true that = 0. If all models fall into this class for pairs with 
Kj_ exactly vanishing, then there should be a wide range of models which are "close" to 
being in this class for pairs with sufficiently small K±. In the next section, we will look at 
an illustrative example model. 



V. A MODEL WITH TRANSVERSE EXPANSION 

Just as longitudinally expanding sources feature z-t correlations in their CM frame, 
transversally expanding sources will in general feature x-t correlations in their CM frame 
or in any frame which is only longitudinally boosted relative to it. Consequently, such 
sources will exhibit nonvanishing B w and/or B 3 i and thus not belong to the class ( [4.191) . 
Nevertheless, by working with Eqs. ( |2.5|) as well as the corresponding general expressions 
for the "GYK" fit, it can be shown that if for some range of K, all of the A M are of the same 
order and 

(AxAq^o) 2 « Pi , (AsA^si) 2 « 0{ , and 1 - « 0[ , (5.1) 
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then Eqs. (|4.21|) are still good approximations. In other words, Eqs. ( |5.1|) define a larger class 



of models for which it is possible to unambiguously extract the parameters Ai(K) ~ X 2 (K), 
As(K), Ao(K), and B 3Q (K). Furthermore, from Eq. Q3.9D it can be verified that if Eqs. ( p.l| ) 
hold, then in the v (K) = frame, AqA^-B^q <C 1, so it is a good approximation to treat 
the curvature tensor as being diagonal, and the correlation radii in this frame are given by 
Eqs. (gig ). 

We will now examine a specific source model which exhibits transverse as well as longitu- 
dinal expansion to see what kind of restrictions Eqs. ( |5.1|) impose on the average momentum 
K. We consider a source function of the form (|4.14j) , but with an expansion four- velocity 
now given by 



u{x) = ( y/l + {v tP /R) 2 diri, (v t x/R), (v t y/R), + {v t p/R) 2 sh V ) . (5.2) 



Since u(x) is well defined and u-u = 1 for arbitrarily large Vt, Eq. (|5.2j) can be used for 
modeling relativistic (v t > 1) as well as nonrelativistic (v t <C 1) transverse expansions. 
(Only for nonrelativistic expansions does the parameter v t represent the transverse velocity 
of the source at p = R.) 

By using Eq. ( |2.2j ) to calculate the saddle point of the emission function and recalling 
that K does not have any component in the y ("side") direction, it is easily found that 
f = Tq and y = 0. Similarly, one can show that 

v t x K\v\jT 



R 1 + (m t v?/T)ch(fj -Y)[l + (v t x / R) 2 ]- 1 / 2 



(5.3) 



and 



! _ M = (v t x/Rnm t v 2 /T)ch(f,-Y) 

A? (m t v 2 /T)ch(fj -Y)[l + (v t x/R) 2 } + [1 + (v t x/R) 2 f 2 ' 

From these equations, it can be seen that if we demand that 



E K <^, (5.5) 

Vt 

then 1 — A|/ Af will always be less than Actually, given the form of Eqs. ( p.3[ ) and 
Q5.4D , we are justified in deducing that Eq. (|5.5|) implies 1 — \\/\\ <C (3j rather than the 
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weaker condition 1 — A^/A^ < 0\. In Appendix B we also show that if (|5.5|) is satisfied, then 
(A1A0-B10) 2 <C 01 and (A3A1-B31) 2 <C 0{. Consequently, in the context of the present model, 
Eq. ( |5.5| ) is a sufficient condition to justify the use of the expressions in the last section. 

In particular, if one makes a "GYK" fit to the correlation function of the model under 
consideration, then for pairs satisfying ( |5.5|) , the extracted correlation radii will measure the 
source in the local longitudinal rest frame at the saddle point. Explicitly, 



Rn 



At 



R t ~ Ai ~ A 2 ^ R 



1 + ^ch(fj - Y) 



T 



-1/2 



Ra 



I^ c h(r7 -Y) - + — 

T Ch[r] Y) ch 2 (f] — Y) (Ar]) 2 



-1/2 



(5.6) 



where r/(K) ~ th _1 [w(K)]. By looking at the m t -dependence of R t , it may be possible to 
extract both the transverse size R and the expansion parameter vf/T. Using the latter, it 
is possible to make a consistency check to see if the pairs under consideration did in fact 
satisfy condition ( |5.5|) . 

The temperature T and transverse velocity parameter vt can also be determined by 
measuring slopes and curvatures of one-particle distributions. In heavy-ion collisions at the 
both the AGS and SPS, these parameters have been estimated to be on the order of T = 140 



MeV and v t = 0.5 pT|-|33|j . From Eq. (|5.5|), these estimates imply that Eqs. ( |5.6| ) should be 
good approximations for pairs with Ek less than about 560 MeV. 



VI. CONCLUSIONS 

We have shown that in general the number of parameters needed to describe a source 
in the quadratic saddle-point approximation exceeds the number of parameters which can 
be determined by making gaussian fits to two-particle correlation data. However, we have 
identified a wide class of interesting models for which the source is fully described by only 
K-dependent parameters which can all be determined from the experimental parameters 
measured in a gaussian fit. Using a realistic three-dimensionally expanding hydrodynamical 
model, we showed that for heavy-ion collisions at the AGS or SPS, it should be simplest 
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to extract the source parameters from correlations of pairs with average energies less than 
about 560 MeV in the measurement frame. It should be noted that the source shape seen 
by these particles may still not be the geometrical shape of the source, but rather that of 
the local region of homogeneity which is affected by the expansion flow profile of the source. 
To separate the flow effect from the underlying geometry, it is necessary to determine the 
K± dependence of the HBT radius parameters. Fortunately, much of this dependence may 
be possible to see while staying in the "simple" regime. For example, for pions with Kl = 
in the measurement frame, values of K± up to 540 MeV/c will still correspond to energies 
below 560 MeV, so at least five 100-MeV bins in K± can be explored below the limit. Of 
course for any given analytic model it should also be possible to extract source parameters 
from large K± correlation radii by using expressions such as those in Eqs. ( |3.3|) . However, at 
large K±, each correlation radius will contain contributions from a large number of effects 
which may be difficult to disentangle. 

In the past there has been some debate as to which longitudinal reference frame would 
be the most appropriate for measuring correlations from a given reaction. For example, for a 
source which is not expanding longitudinally, the source center-of-mass frame is the natural 
choice. On the other hand, for an infinite source which is undergoing a boost-invariant 
expansion, the LCMS represents the local rest frame of the source and is thus the natural 
choice. Since experimental reactions undoubtedly produce sources which lie somewhere 
between these two extremes, some intermediate frame is needed. The generalized Yano- 
Koonin fitting procedure of Eq. (|3.7|) has the advantage that it does not require one to 
postulate a reference frame beforehand; the data themselves determine a frame for each 
value of K. For many interesting "intermediate" models, the parameters in this v(K) = 
frame measure the source in its local longitudinal rest frame, while fixed frame or LCMS 
parameters measure the source in some different frame. 

Certainly the best way to compare any given model to correlation data is to make the 
comparison directly in a six- dimensional (pi,P2) space, rather than to compare the fitted 
correlation radii of the model to those extracted from the data. Nevertheless, we have shown 
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here that gaussian fits can still reveal some very interesting information about the velocity 
of the local longitudinal rest frame as well as the lengths of homogeneity of the source. 
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APPENDIX A: 

In this appendix we show how the $ dependencies in Ao and B 30 of Eqs. (|4.21|) can be 
removed through the introduction of new fitting procedures. We label the fitting procedure 
defined by Eq. (|3.2| ) with an "a". Thus, R 2 (a) re fe rs to the square of the "longitudinal" 
radius as found by fitting the correlation with Eq. ( |3.2| ). We define a "6" fitting procedure 
by using (3iq 3 instead of 53 in making a gaussian fit to the correlation function. In other 
words, 

C(q,K) = 1 ±Xex V [-q 2 1 R 2 1 (b) - q 2 2 R 2 2 (b) - (fcq,) 2 R 2 3 (b) - 2g 1 (/3 1 g 3 )#? 3 (&)] , (Al) 

where we have suppressed the K dependence of the R 2 parameters. It can now be seen that 
in the LCMS (primed) frame 



[iff (q) - iff (a)] iff ( q) - [iff 3 (a)f 
iff(6) 

[Rf(a)-Ri(a)]R'i(a)-[R? 3 (a)f 



A, 
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B' = , . (A2) 



To calculate Ao and -B30 in a fixed frame, we need to introduce two additional fitting 
procedures. Procedure "c" is defined by using foqi, P 3 q2, and /3iq 3 in place of q%, q%, and 
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g 3 when making gaussian fits, while procedure "d" is defined by using J\(3z\q\, y | /?3 1 ^2 , and 
qz- We then have 

2 _ [R 2 (a) - R 2 (a)} R 2 (a) - [R 2 3 (a)} 2 

Rl(c)-2RUc)+Rl(c)-Ri(c) 
B RUb)TjRl(d)-R 2 2 (d)] 

30 [R\{a)-Rl{a)]Rl{a)-[RUa)?' 
where the — (+) sign in B 30 refers to bins in which /3 3 > ({3$ < 0). Note that due to this 
distinction, the above method for determining source parameters should not be applied to 
bins in which some of the pairs have (3% > while others have (3% < 0. 



APPENDIX B: 



To prove the remaining inequalities in Eq. ( |5.1|) , we begin by presenting the following 
easily verified relations: 



B 



m t v t (v t x\ ( t\ sh(?? - Y) 



RT \R ) \r 2 ) h + (v t x/Ry 



TTltVt (v t x\ ( z \ sh(i] — Y) 
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Bl ° RT V R) \tI) ^1 + (v t x/Rf 
1 m tv 2 t ch(r/-F) 

\\ R 2 T [1 + {vtx/Rff/ 2 ' 1 ; 



The saddle point in spacetime rapidity for this source is given by 



th(fj - Y) - y/1 + {v t x/RY sh(77 - y) _ = . (B2) 

To derive the needed expressions for Aq and A§, we have first used Eq. ( P2|) to check numer- 



ically that if {m t /T)Jl + (v t x/R) 2 > 0.7, then for F ^ 0, \r]\ < \Y\ and r\ has the same 



sign as Y. Although (m t /T)yl + (v t x/R) 2 > 0.7 may present a significant restriction for 
electron or photon correlation measurements, it does not present a significant restriction for 
current two-hadron correlation measurements from particle or heavy-ion collisions. Even for 
pions, temperatures of up to 200 MeV would still satisfy this condition. Given this condition, 
from Eq. ( p2|) it is possible to show that for Y ^ 0, 
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Using this inequality, then for all Y it is true that 

2 „1,2, 



The desired inequalities are now easily proven: 

(AxAa^i) 2 < (v t x/R) 2 < & 2 
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